Theory for Magnetism and Triplet Superconductivity in LiFeAs 
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Superconducting pnictides are widely found to feature spin-singlet pairing in the vicinity of an an- 
tiferromagnetic phase, for which nesting between electron and hole Fermi surfaces is crucial. LiFeAs 
differs from the other pnictides by (i) poor nesting properties and (ii) unusually shallow hole pock- 
ets. Investigating magnetic and pairing instabilities in an electronic model that incorporates these 
differences, we find antiferromagnetic order to be absent. Instead we observe almost ferromagnetic 
fluctuations which drive an instability toward spin-triplet p-wave superconductivity. 



PACS numbers: 74.70.Xa, 74.20.Rp, 74.20.Mn 

I. INTRODUCTION 

The family of pnictide superconductors has generated 
large interest during the last two years, because their 
superconductivity (SC) is generally thought to be un- 
conventional, i.e., not due to phonons, and because some 
pnictides have high transition temperatures reminiscent 
of cuprates. The details of the pairing are far from being 
clarified, and may vary significantly between the different 
families; for two recent reviews see Refs. 11|2 , A plausi- 
ble scenario for the pairing mechanism involves exchange 
of virtual spin fluctuations, as SC is often observed in 
the vicinity of an antiferromagnetic (AF) spin-density- 
wave (SDW) state. In the 1111 and 122 families, the 
SDW is apparently driven by excellent nesting of electron 
and hole Fermi pockets. When doping reduces the nest- 
ing, SDW order is weakened, but AF fluctuations remain 
and are believed to provide the pairing for the Cooper 
pairsPS Together with the high transition temperatures 
observed in some pnictides, such a scenario raises hopes 
that pnictides might shed light onto similarities - or dif- 
ferences - between spin-fluctuation superconductors like 
heavy-fermion systems, with typically lower transition 
temperaturespEl and the more enigmatic cuprates. 

The recently-discovered pristine superconductor 
LiFeAs differs in key respects from the most commonly 
studied pnictide families. Doping is not necessary for 
superconductivityP^the nesting between hole and elec- 
tron pockets is rather poor and consequently an SDW 
is not observed. Moreover, compared to other pnictide 
families, LiFeAs has much shallower hole pockets around 
the center of the Brillouin zonePHXI It was pointed out 
that the flat top of these pockets implies a large density 
of states,^ which one can expect to promote ferromag- 
netic (FM) fluctuations. Indeed, recent nuclear magnetic 
resonance (NMR) measurements indicate a relatively 
high susceptibility within the a-b planes. Remarkably, 
the Knight shift does not change at the superconducting 
transition for magnetic fields perpendicular to the c 
axisp^l which is evidence for triplet superconductivity. 
This discovery relates a member of the pnictide family 
to yet another family with unconventional SC, the 



ruthenatesP^ 

In this paper we model the distinguishing electronic 
features of LiFeAs, i.e., the absence of nesting and the 
presence of shallow hole pockets. We examine the mag- 
netic properties of this Hamiltonian by combining a 
mean-field approximation and the random-phase approx- 
imation (RPA). We find that the magnetic fluctuations 
are dominated by almost FM processes across the smaller 
hole pocket. We proceed to investigate what types of 
Cooper pairs can be stabilized by the spin and charge 
fluctuations. We use the same approach that has been 
used to treat other pnictides, and which yields AF order 
and singlet pairing in the 1111 and 122 familiesP^^ In 
the case of LiFeAs, however, we find that the almost FM 
fluctuations can drive triplet p-wave superconductivity. 
Throughout, we discuss the connection of our work to 
experimental results. 

The possibility of triplet pairing was already discussed 
in the early days of pnictide research. The degeneracy of 
the xz and yz orbitals that dominate the Fermi surface 
(FS) permits spin-triplet and orbital-singlet pairs with 
(nodal) s- or d-wave-like character QSH221 Triplets with 
p-wave symmetry, stabilized by large Hund's rule cou- 
pling, were reported for a three-band model and involved 
mostly the electron pocketsP^l The mechanism proposed 
in this paper, however, is very different and specific to 
LiFeAs, originating from the high density of states due 
to the rather flat 'roofs' of the shallow hole pockets. 

II. MODEL 

In order to capture the most relevant states around 
the FS, we model the two concentric hole pockets and 
the two electron pockets, which can be achieved using 
three orbitals per iron ion. We use a one-iron unit cell, 
with k running over — tt < k x , k y < tt in the correspond- 
ing extended Brillouin zone, and the hole pockets are 
then around T = (0,0), while the electron pockets are 
found at X — (tt, 0) and Y — (0,tt). The three orbitals 
with the largest weight in these three pockets are the 
xz, yz, and xy orbitals, denoted by indices 1, 2 and 3. 
The momentum-dependent non-interacting Hamiltonian 
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FIG. 1: (Color online) (a) Band structure of (dashed lines) the 
three-band model Eq. |TJ from Ref. and (solid lines) the 
model for LiFeAs. (b) Fermi surface and dominant orbitals for 
the LiFeAs model; ai and a.2 are hole pockets, /3i and 02 are 
electron pockets. Parameters: ti = 0.02, ti = 0.12, tz = 0.02, 
U = -0.046, t B = 0.2, t 6 = 0.3, t 7 = -0.15, t 8 = -t 7 /2, 



t 9 



-0.06, t 10 



-0.03, tn = 0.014, A,, 



1, fx = 0.338. 



The filling is four electrons per site, energies are in eV. 



is given by 



(1) 



where „ a (c? k CT ) creates (annihilates) an electron in 
orbital /i with momentum k and spin a. The elements of 
the hopping matrix ^'"(k) are 

yii/22 _ 2t 2 /\ cos k x + 2tx/ 2 cos k y + 4i 3 cos fc a cos k y 

±2tn(cos2k x ~ cos2k y ) ~ [i, (2) 
T 33 = A xy + 2t 5 (cos k x + cos k y ) + 4i 6 cos k x cos fc y 
+ 2tg (cos 2fca; + cos 2k y ) 

+ 4fi (cos 2k x cos k y + cos fc^ cos 2k y ) — fx, (3) 

T 12 = ^21 = _ 4<4 sin fc ^ gin ( 4 ) 

T 13 = T 31 = 2it 7 sin k x + 4it 8 sin k x cos , (5) 
T 23 = T 32 = 2it 7 sin k y + Ut % sin fc y cos k x , (6) 

where a bar on top of a matrix element denotes the com- 
plex conjugate. The parameters are given in the caption 
of Fig. [TJ 

The one-particle bands defined by this kinetic energy 
are shown in Fig. 1(a) for two parameterizations: one is 
the three-band model of Ref. [351 where only hopping to 
first and second neighbors is included. While this model 
already has rather shallow hole pockets compared to ab 
initio bands for most pnictide compounds, it can repro- 
duce spectral feat ures o f the AF phase observed in many 
undoped pnictide j2512SI and has an instability to singlet 
pairing.U-3 The second more refined model makes use of 
longer-range hoppings in order to describe the bands of 
LiFeAs more closely. The two hole pockets, which have 
mostly xz and yz character and are degenerate at T, are 
far shallower than the electron pockets and have very 
different radii in accordance with ARPES findings and 
ab initio density functional calculations 12HJJJ Both mod- 
els reflect the poor nest ing of the electron and hole pock- 
ets in LiFeAs, see Fig. |l(b) As will be discussed later, 
the flat top of the hole pockets has a crucial impact 
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FIG. 2: (Color online) Magnetic phase diagram of the three- 
band model for LiFeAs. Inset: Magnetic susceptibility x zz (q) 
for U = 0.42eV and T = 35 K and T = 135 K. 



on the model properties; angle-resolved photo-emission 
(ARPES) dataP might suggest even flatter hole pockets. 

We include all Coulomb interactions between electrons 
in the same Fe atom: 

Hi = U^rii^f n i>Vti , + V ^ X! n 'w n i,n,<r' ( 7 ) 

\.v ijU^^i er, a' 

Here n\^ v _ a (S^,,) is the number (spin) operator for the 
orbital v on site i. Invariance under rotation of the orbital 
degrees of freedom is ensured by setting J = J' and V = 
(2U - 5J)/4. We only present results for J = 0.25U, but 
we have verified that other choices lead to similar results. 



III. MAGNETISM 

Following Ref. [STJ we construct the magnetic phase 
diagram by first using a mean-field ansatz restricted to 
FM and two-site AF phases^. Unlike models of the 1111 
and 122 pnictides, here only FM and PM phases are sta- 
bleP3 The phase diagram is then refined by evaluating 
the static transverse spin susceptibility x (q.) within 
the RPA. This yields two important results: firstly, nega- 
tive values of x (q) signal an instability toward a mag- 
netically ordered state (FM for q = 0, AF otherwise); 
secondly, in the PM phase the maximum of x ^ (q) indi- 
cates the wavevector of the dominant spin fluctuations. 
In the PM phase x h ( f l) — 2x zz (q) is obtained by solv- 
ing the Dyson equation x S = + X U S X S for the 
irreducible spin susceptibility x S29 . The non-zero ma- 
trix elements of U s can be found in Refs. 14 16, The 
Lindhard function x^ is defined 



3) / . \ 1 \ - \ - n F (E jM ) - n F (E fM+q _) 

x u J> ,(k)^ >(tt (k)u i , iAl ,(k + q)M^ >I/ (k + q) (8) 
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where v, 1/ , fi, // refer to the orbital, np(E) is the 
Fermi function, and Ej,± are the eigenvalues of HqQs.) 
[Eq. (IT l]. The coefficients Uj „(k) transform the diagonal- 
izing (band) annihilation operators 7. k (corresponding 
to Bj,k) into the orbital basis, i.e. cZ„ k = J2j u j,vQ*)lj k- 
The static transverse spin susceptibility is then written 

x~ + (q) = £„, /1 xfU Ml/1 (q.w = °)- 

Figure [2] shows the phase diagram as a function of U 
and temperature T. At low t/, the model remains PM 
down to T = K, but the dominant fluctuations change 
from FM to incommensurate (IC) at T w 100 K. As [/ 
is increased, the IC fluctuations can drive a transition to 
an IC-AF state, while at higher U a FM phase is real- 
ized. The wave vectors characterizing the IC fluctuations 
in the low-T PM state lie on a ring-shaped feature cen- 
tered at q = (0, 0), as can be seen in the example plotted 
as the inset in Fig. [2] The radius of the ring is exactly 
twice the radius of the inner hole pocket, revealing that 
it originates from scattering diagonally across this FSpSS. 
This scattering is strongly favored because of the identi- 
cal orbital composition of diagonally separated parts of 
theFS. 

Experimentally, a relatively large susceptibility is ob- 
served parallel to the a6-plane, which slightly decreases 
with T. 12 While an increasing susceptibility is expected 
for dominant FM fluctuations, this observation is consis- 
tent with 'almost FM' IC fluctuations, where the weight 
observed at (0,0) in x zz ( ( l) no longer grows with lower- 
ing T below the transition from FM to IC fluctuations. 
It may even slightly decrease as the ring feature becomes 
more pronounced, and our model shows such behavior at 
low temperature T < 50 K. These results indicate that 
LiFeAs is relatively near to a FM instability, which can 
be triggered by reducing the size of the inner hole pocket, 
e.g. by electron doping. 



IV. PAIRING SYMMETRIES 

The strong spin fluctuations in the low-T PM phase 
may drive the pairing of electrons. To determine the 
symmetry of a possible superconducting state, we employ 
a weak-coupling method due to Scalapino et alW^ which 
has been widely used to study the pnictides!^^ 

The pairing vertex due to the exchange of spin and 
charge fluctuations is obtained within the fluctuation ex- 
change approximation. We have 
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FIG. 3: (Color online) Triplet pairing, (a) Eigenvalues A 
of the leading gap symmetries as a function of U. (b) Gap 
around the Fermi pockets for the p x -w&ve pairing, which is a 
member of the degenerate pair that gives the dominant triplet 
p-channel. The winding angle about the FSs 9 is measured 
with respect to the fc^-axis, taken in the anti-clockwise di- 
rection, (c) Gap magnitude of pairing states combining the 
degenerate p x and p y solutions of Eq. l|12| ). ai, oli, /3i, and 
02 refer to the FS pockets, see Fig. 1(b) 



for singlet and triplet pairing, respectively. Here \ c is 
the RPA irreducible charge susceptibility, which obeys 



the Dyson equation x C = — X^U c x C ■ The non- 
zero elements of U are again found in Refs. IMHTBI 

Assuming that the dominant scattering occurs close to 
the FSa^EII we describe the scattering of a Cooper pair 
from the state (k, — k) on FS Cj to the state (k', — k') on 
FS Cj by the projected vertices 

r£7 s '*(k,k') = ^7(-kH,«(kK^(k>* 5 (-k') 

a, /3, 7,<5 

xRe{r^-(k-k', W = 0)} . (11) 

The superconducting gap on the FSs is written as A(k) = 
Ag u (k), where g v {}s) is a dimcnsionless function describ- 
ing the momentum dependence, and v = s {€) denotes 
a singlet (triplet) state with even (odd) parity. g v (\s) is 
obtained by solving the eigenvalue problem 

r dk' 

- E 4 4^ko r - (k ' k '^ (k,) - x ^ k) (12) 

where vp^(k) is the Fermi velocity and A is a dimension- 
less coupling strength. The gap function g v {k) with the 
largest A has the highest T c , and hence fixes the symme- 
try of the pairing state. We use a 192 x 192 k-point mesh 
and T — 50 K to determine the leading pairing instabil- 
ity. Working at T < 50 K requires a much larger k-point 
mesh, especially when small momenta are important for 
the pairing. 



Figure 3(a) shows the eigenvalues A of Eq. (12 1 as a 
function of U for various pairing symmetries at T — 50 K. 
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The dominant pairing channel corresponds to a triplet 
gap with p-wave symmetry. The gap function g x (k) 
for the p x -wave state is shown in Fig. |3(b)| this state 
is degenerate with a pj,-wave state with gap function 
9y{k Xl k v ) = g x {k y ,k x ). In order to maximize the gap 
at the FS, these two states can be combined to form a 
number of unitary states: (i) d(k) = A[g x (k)±ig y (k.))e z , 
(ii) d(k) = A[g x (k)e x ± g y (k)e y ], and (hi) d(k) = 
A[g x (\t)e y ± g y (k)e x ]. All of these states have the same 
gap magnitude, shown in Figure 3(c) and are hence 
degenerate in the present model. The inclusion of ad- 
ditional interactions, e.g. spin-orbit coupling, ma y lift 
this degeneracy and favor one state, as in S^RuOj 13 * 32 ! 
The absence of the Knight shift for fields parallel to the 
ab plane, and the strong out-of-plane anisotropy, sug- 
gests that the opposite-spin-pairing d(k) = A[g r a; (k) ± 
ig y (h)]e z state is most likely in LiFeAsP^l 

The sub-dominant pairing states are a triplet /-wave 
state and a singlet I — 6 state. The latter can become 
rather strong close to the IC-AF phase, where the small-q 
IC fluctuations favor a singlet gap that changes sign over 
the same small momentum difference. It is interesting to 
compare our results with the related three-orbital model 
proposed in Ref. 1251 where scattering across the inner 
hole pocket also produces a similar ring-shaped feature in 
X zz (q) ■ The radius of this hole pocket is much larger than 
in our model for LiFeAs, however, and there is hence not 
enough weight at q = (0, 0) to support triplet pairing.^ 
Instead, the larger radius of the ring feature favors slowly- 
varying s-wave or B\ g singlet gaps. Since triplet SC is 
dominant as long as the weight at q = (0, 0) remains high 



enough down to the critical temperature, we conclude 
that it will be favored by a small inner hole pocket. Other 
features of the bands play only a secondary role, e.g. 
by slightly varying the hopping integrals we can find a 
much smaller gap on the inner hole pocket, or an almost 
constant gap on the electron pockets apart from very 
narrow nodal regions. 



V. CONCLUSIONS 

In this paper we have presented a coherent picture 
for the properties of LiFeAs driven by magnetic fluctu- 
ations. Due to the poor electron-hole nesting and the 
shallow hole pockets, we find that scattering across the 
small inner hole pocket dominates the physics, leading 
to a ring-shaped feature around (0, 0) in the momentum- 
dependent magnetic susceptibility. Due to these almost 
FM processes, triplet p-wave superconductivity is favored 
over singlet pairing. 
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